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Abstract 

We calculate axion radiation emitted in the collision of two straight global strings. The strings 
are supposed to be in the unexcited ground state, to be inclined with respect to each other, and 
to move in parallel planes. Radiation arises when the point of minimal separation between the 
strings moves faster than light. This effect exhibits a typical Cerenkov nature. Surprisingly, it 
allows an alternative interpretation as bremsstrahlung under a collision of point charges in 2 + 1 
electrodynamics. This can be demonstrated by suitable world-sheet reparameterizations and di- 
mensional reduction. Cosmological estimates show that our mechanism generates axion production 
comparable with that from the oscillating string loops and may lead to further restrictions on the 
axion window. 

PACS numbers: 11.27.+d, 98.80.Cq, 98.80.-k, 95.30.Sf 
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I. INTRODUCTION 

For more than twenty years the Peccei-Quinn axion remains one of the most 

serious dark matter candidates. Its properties depend essentially on one unknown parameter, 
the vacuum expectation value / marking the energy scale of the U(l) symmetry breaking. 
The axion acquires a small mass 

lO^^GeV 



6 ^eV 



f 



(1) 



in the QCD phase transition. The quantity / is a free parameter of the theory, but it is 
severely constrained both by the accelerator physics and the astrophysics, the upper bound 
from the combined data being several meV. The lower bound is of cosmological origin and 
follows from the requirement that axions produced during the cosmological evolution do 
not overdose the Universe Q]. The cosmological axion production is mainly due to the 
axion radiation from global strings and domain walls. In view of still persisting theoretical 
uncertainty, the corresponding lower bound lies between several units and several tens of 
fieV, corresponding to the maximal value of / between 10^^ and 10^^ GeV. (For recent review 
of the present theoretical status of the axion see |5|; the current astrophysical status of this 



model is reviewed in 



earlier reviews include 



The axion string network IJ, ll^ ll^ is formed at the temperature of 

the Peccei-Quinn phase transition /, and it is usually assumed that the reheating tempera- 
ture is higher than /, otherwise the network would get diluted by the exponential expansion; 
this introduces more uncertainty into theoretical predictions. Strings are primarily produced 
as long straight segments with length of the order of the horizon size, and they initially move 



with substantial friction 



0,0, 



19| due to scattering on them of the cosmic plasma particles. 



At some temperature T* < f the scattering becomes negligible and the string network enters 



the scale-invariant regime 
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when strings form closed loops via interconnections 



23 



. The standard estimates of the 



and move almost freely with relativistic velocities 
axion radiation from global strings are based on the assumption that the main contribution 



comes from the oscillating string loops 
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28|, 



29 



Sd, ^ The amount of 



radiation is large enough and leaves a rather narrow window for possible values of /. 

Here we would like to discuss one additional mechanism of axion radiation: the 
bremsstrahlung which can be produced in collisions of long strings. In fact, the number 
of collisions of long strings is at least not smaller than the number of string loops which are 



formed under some of these collisions (note that in our case passing of strings through each 
other is not assumed). Meanwhile, as far as we are aware, this mechanism was not explored 
so far in the axion physics context. It works as follows. Neglecting the boundary effects, 
let us assume that two infinite straight strings move in parallel planes, the two strings be- 
ing generically inclined with respect to each other. Then the point of minimal separation 
between the strings, which marks the localization of the effective source of axion radiation, 
is allowed to move with any velocity, superluminal in particular. In this latter case one can 
expect Cerenkov axion emission to emerge. Similar mechanism was earlier suggested for 
gravitational radiation of local strings [2^, but in that case the explicit calculations have 
led to zero result. Vanishing of the gravitational radiation, however, has a specific origin 
related to the absence of gravitons in the 2 + 1 gravity theory. Indeed, as it was explained 
in js^, two crossed superluminal strings can be "parallelized" by suitable coordinate trans- 
formations and world sheet reparameterizations, so the problem is essentially equivalent to 
the point particles' collision in 2 -j- 1 gravity theory. For other fields, (e.g. electromagnetic) 
this mechanism does work; see [s^l for calculation of the Cerenkov electromagnetic radiation 
generated in collisions of superconducting strings. So one can expect to have a non-vanishing 
axion bremsstrahlung from collisions of global strings as well. 

Our calculation follows the approach of 3^ and it is essentially perturbative in terms of 
the string-axion field interaction constant (equal to /) involving two subsequent iterations in 
the strings equations of motion and the axion field equations. This is similar to calculation of 
the bremsstrahlung from charged ultrarelativistic particles. In this latter case the iteration 
sequence converges if the scattering angle is small. We assume the same condition to hold 
for the strings, though we do not enter here into a detailed investigation of the convergence 
problem. We would like to mention also a similar perturbative approach in terms of geodesic 



deviations 



3- 



II. GENERAL SETTING 

Consider a pair of relativistic strings 

x'^ = <(0, /x = 0,l,2,3, aa = {T,a), a = 0, 1, 
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where n = 1,2 is the index labelhng the two strings. The 4- dimensional space-time is 
assumed to be flat and the signature is chosen as +, and (+, — ) for the string world- 
sheets. Strings interact via the axion field as described by the action [l^ 



S = Sst + SB, (2) 



where the string term is 



^^t = -Yl [ {^V^r'dax':,d,xlr,^, + 27r/„S^,e'^^9,x^:9fex::)rfV„, (3) 

^ — 1 o ^ 



Sb = \ I H.^xH^'^^d^x. (4) 



n=l,2 

and the field action is 

6 

Here are the (bare) string tension parameters, /„ are the corresponding coupling parame- 
ters (their n-labelling helps to control the perturbation expansions), the Levi-Civita symbol 
is chosen as = 1, 'jab is the metric on the world-sheets. 
The totally antisymmetric axion field strength is defined as 

H^u\ = df^B^x + duB\^ + d\B^y. (5) 

Variation of the action Q over x^ leads to the equations of motion for strings 

lld^iV^l'^'d.xO = 2vr/„if^^C^ (6) 

where 

Vf = e'^'daxy.xt (7) 

and H'^p is the total field strength due to both strings (no external field is assumed). The 
corresponding potential two-form is the sum 

Bi"" = 5^ + B^\ (8) 

of contributions B^'-' due to each string: 

= ^rrXar (9) 

Here the Lorentz gauge d\B^^ = is assumed, the 4-dimensional D'Alembert operator is 
introduced as □ = —rj'^^d^dy, and the source term reads 



The self-action terms in the equations of motion (jH)) diverge both near the strings and at 
large distances, so two regularization parameters 6 and A have to be introduced |37,] . These 
can be absorbed by classical renormalization of the string tension as follows 

fi = fio + 27if\og{A/6). (11) 

The ultraviolet cutoff length 6 is of the order of string thickness 5 ~ 1//, while the infrared 
cutoff distance A can be taken as the correlation length in the string network. 

Assuming that such a renormalization is performed, we are left with the same equations of 
motion (jH)) with the physical tension parameters yU„ which contain only the mwtwa/ interaction 
terms on the right hand side. The constraint equations for each string read : 

daX^dhX-'V.u - \labl'\dcX^d<lX''^^.u) = 0. (12) 

Using the invariance of the action (jSJ under separate world-sheet reparameterizations 

< ^ KK) (13) 

and under Weyl rescaling of the metric 'jat, one can fix the flat gauge 

lab = Vab, (14) 

SO that the constraint equations simplify to 

i;2 + x'^ = 0, xf^x'^r]^^ = (15) 

for each string, where dots and primes denote the derivatives over r and a as usual. In the 
gauge (jl4j) the renormalized equations of motion read : 

/i„ (9,2 - = 27r/„/7„TXA, n^n', (16) 

where the field H^^^ corresponds to the contribution of the n'— th string (no sum over n'). It 
is worth noting that imposing flat world-sheet metric we still do not fix the gauge completely, 
since one is free to perform two-dimensional linear transformations preserving the Minkowski 
metric ()14|) on the world sheets. 

The retarded solution to the wave equation can be presented in its standard form 

B^" = ^Ti j G'^\x - x')J^%x')d^x' , (17) 



where the Green's function is 

G'^'^x - x') = / 4 (18) 

with e —>■ +0. The advanced Green's function is obtained by changing the sign of e. 
The energy momentum tensor for the axion field reads 



(19) 



To avoid comphcations with definition of the wave zone for infinite strings, it is convenient 
to compute the radiation reaction power using the divergence equation 



-AnH^'^^J^p, (20) 



dx^ 

where H^""^ should be taken as the half- difference of the retarded and advanced solution of 
the wave equation, H'^.^f = (i^^dt^ ~ ^bA^)/'^- '^^^ corresponding potential has then the 
following Fourier transform: 

W\ 



B^^^{k) = -At7r' — J^,(k)5{e) , (21) 



where the current 



is transverse 



J^Ak) = I J^,uix)e"'^d''x (22) 



k^J^.u{k) = (23) 

and taken on the mass-shell A;^ = of the massless axion. The radiative loss of the 4- 
momentum during the collision can be evaluated as 

= j d^Tf^'d^x = -An J HZ^J^pd^x . (24) 

Substituting (jHl) into (j21I) we obtain 

= II k''^^\Jap{k)\'S{e)d'k. (25) 

Alternatively, this quantity can be presented as a square of the polarization projection of 
the Fourier-transform of the current. Indeed, in three dimensions the axion field propagating 
along the wave vector k, k^ = {u, k), has a unique polarization state 

^ e.^e^-efe?), z,j = 1,2,3, (26) 



where and e'^ are two unit vectors orthogonal to k and to each other: 

e"^ ■ = 0, k ■ e, = 0. (27) 
Using antisymmetry and transversahty of the current (j^!^ and the completeness condition 

e'^e'^+e'te'^ = 6,,-k,k,/u', (28) 

one finds 

AP^' = - I k^^AX^{k)eiA^6{k^)(fk. (29) 
Integrating over /c°, we finally obtain 

1 C kf^ 

AP^ = - / ^|J(k)|2d3A;, (30) 

TT J |k| 

where 

J(k) = r^{k)e,j (31) 

with = |k|. In what follows we shall use the following explicit parameterization of three 
orthogonal vectors by spherical angles : 

k = co'[sin6' cosy?, sin^sinyj, cos^], 

e"^ = [— siny9, cosy^, 0], (32) 
= [cos cosy?, cos 6' sin yj, —sin 6*]. 

Our strategy consists in solving the system of equations (fTBl ITTj) iteratively, by expanding 
all functions in powers of the coupling constant /: 



oo 



= (33) 

1=1 

CO 

J^'^ = ^iJ^", (34) 

1=0 

oo 

: $^/x^(a). (35) 



1=0 

Such expansions have to be performed separately for each string. In our notation, the 
zeroth-order approximations for the string embedding functions Qx!^{a) give zeroth-order 
source terms qJ^'^ , which generate the first order fields iB^'^ originating on each string. 
These quantities, when substituted to the string equations of motion (fTH|l. will generate the 



FIG. 1: In the chosen frame, the first string is at rest, the second moves in the parallel plane 
separated from the first by an impact distance d = d2 — di being inclined at the angle a with 
respect to the string at rest. Its velocity v is perpendicular to the string itself. 

first order terms ^^x^, and so forth. Convergence of this expansion (with a formal parameter 
/ which has a dimension of an inverse length) is difficult to explore in general. This is likely 
to depend substantially on the choice of the zeroth-order solution of the string equations of 
motion. In this paper, the strings are assumed to be in the unexcited state in the zeroth 
order, that is they are straight and freely moving. In this case it can be expected that 
the power series solutions will give sensible results (at least in the lowest appropriate order 
when radiation appears) if the deviation angle of colliding strings due to their lowest order 
interaction is small. We shall come back to this condition later on. 

III. COLLISION KINEMATICS 

The kinematics of collision is shown on Fig. 1. The unperturbed world-sheets of un- 
excited, freely moving straight (infinite) strings are described by linear functions of r, a 

ox^ = < + <7:„ +S;^a„, n = l,2, (36) 

with c?^ and being space-like and time-like constant vectors. Then the constraint 
equations imply the orthogonality and normalization conditions (choosing the unit norm 
time-like 4-velocity ) 

M^S^ = 0, M% = l = -WY.^ (37) 



8 



for each string. We choose the reference frame such that the first string is at rest and is 
stretched along the axis: 

<= [1,0,0,0], S5' = [0,0,0,1]. (38) 

The second string is assumed to move in the plane x"^, x^ with the velocity v orthogonal to 
the string itself: 

^2 = 7 [1, 0, — f cosa, f sina]. Eg = [0, 0, sin a, cos a], (39) 

where 7 = (1 — f ^)^^/^. We also choose both impact parameters di^ to be orthogonal to 
and and aligned with the axis x^, the distance between the planes being d = d2 — di. 
The angle of inclination a of the second string with respect to the first one can be presented 
in the Lorentz-invariant form as 

a = arccos(— S1S2), (40) 

and the relative velocity of the strings as 

V = {l-{uiU2)-^)'^. (41) 

With this parametrization of the unperturbed world-sheets, the point of minimal separation 
between the strings (which always remains equal to d) moves with the velocity 



_ — ^ ; I (42) 

sma \1 - [2:12:2)^ J 

along the x^-axis. This motion is not associated with propagation of signal of any kind, so 
the velocity Vp may be arbitrary, in particular, superluminal. The case of parallel strings 
corresponds to = 00. 

Let us explore the residual gauge invariance consisting of the S0(1,1) transformations of 
the world-sheets 

Tn r^coshxn + c^nsinhxn, 

cr„ r„sinhxn + acoshxn, (43) 
preserving the constraints p5|) . It is easy to show that the relative orientation of two strings 



)2-l\^ 



is not a gauge independent property 



One can build the following matrix describing the 



relative orientation of the world-sheets 

Kab = dax'^dbx'^rinu, (44) 



whose determinant 

K = det Kab = 7 cos a (45) 

is invariant under the transformation (021) • But the relative velocity (PT|) and the inclination 
angle (|4U|) separately are not invariant. In particular, for superluminal strings one can always 
find the world-sheet boosts which make the strings parallel. This is achieved by the following 
choice of parameters in ()43|) for two strings 

sin a tana 
tanhxi = , tanhx2 = , (46) 

V jv 

which exists if fp > 1. Note that strings are superluminal in the above sense provided that 

K = 7 cos a > 1. (47) 
If this condition is not fulfilled, then the "parallelizing" transformation does not exist. 

IV. FIRST ORDER INTERACTION 

According to the conventions chosen above, the field equations (jHI) for the first order field 
terms ^B^'^ read : 

□ ,Sr = 47roJr, (48) 



with the source terms (14? 

o^r = Y / VrS' (x - oX„(a„)) dV„, (49) 
where the zeroth-order quantities are found under the right-hand side integral : 

vr = <K - K<- (50) 

In what follows we pass to the momentum representation according to the convention : 

F„)^/fWe-.V (51) 

where qx = q^x^. We will use similar notation for all scalar products, when there is no risk 
of ambiguity. The Fourier transform of the zeroth-order current then reads : 

oJr(g) = 27rVne^'''"5(gUn)5(gS„)Kr- (52) 
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The retarded solution iB'^'^ of the Eq. generated by the zero-order source is 



" g2 + 2zeg0 ' ^^^^ 



Note that two 5-functions in this expression shift the momentum q from the pole = 0, so, 
as can be easily seen from the Eq. (jS^I), the first order field is non-radiative. 

The next step is to find corrections to the world-sheet embedding functions due to inter- 
action between the strings via the axion field. Substituting ()53p into the right hand side of 
the equations of motion ([1^) we obtain the following first-order perturbations of the string 
world-sheets 



where 



-Bi = (gMi)(M2Si) - (Sig)(MiM2), Ci = (MiS2)(Sig) - (gMi)(SiS2) . (55) 



Similarly for the second string 

xi^ = I / Xi; — d q, (56) 

J q^[{q^2r - (gM2)2l 



where 



= q^A + B2T.t + C2<, B2 = {qu2){uj:2) - (S2g)(MiM2), 

C2 = (M2Si)(S2g)-(gM2)(SiS2), (57) 

with the same A. It can be checked that the quantities X!^ satisfy the conditions 

X^ui^ = X^Si^ = X^U2^ = X2 S2^t = 0, (58) 

which ensure the fulfilment of the constraint equations (fTH|) up to the first order. 

To estimate the scattering angle 5ip of the string collision let us consider the case of 
parallel strings, a = 0. In the chosen frame one can use as an estimate for 5(p the ratio of 
the limiting value of the component of the 4- velocity of the second string and the initial 



x'^ component: 



6^ ~ lim i^. (59) 
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To calculate the quantity in the numerator one has to integrate in ()56p over d'^q. Integration 
over g° and is performed using delta-functions. Then, in the limit r ^ oo, an asymptotic 
relation 

hm = ,S{qu,) (60) 

T^oo qU2 

gives one more delta- function, allowing integration over g^. The last integration over q^ 
reduces to 



(61) 



the resulting estimate for the scattering angle being 



- (62) 
where we have set /i = /2- Using the Eq. (1111) f or the renormalized string tension, in which 



we assume the second term to be dominant 



221, we obtain 



log(A/d)7f'' 

The smallness of this quantity may serve as a rough estimate of the validity of the per- 
turbation expansions Therefore our iteration procedure works in the relativistic 
case f ~ 1, especially for ultrarelativistic collisions, 7^1. But for global cosmo logical 
strings the quantity log (A/5) is large enough. This improves the convergence of our iter- 
ation scheme even for 7 ~ 1. Therefore our approach seems to be reasonably justified in 
reahstic cosmological situations. 



V. RADIATION 



Using the first order corrections to the world-sheets of the string due to their lowest 
order interaction terms one can evaluate higher order terms of the expansion. Radiation 
emerges in the second approximation for the axion field generated by the source term 
iJ^'^ which will be obtained from the first order contributions to the world-sheet embedding 
functions as follows: 

(5^(x - oa;„((T„)), (64) 



nr-, 1 O ^ 



n 
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where brackets denote anti-symmetrization over indices with the factor 1/2. This follows 
from the fact that the Fourier-transform of this current is non-zero on the axion mass-shell 

The Fourier transform of the first-order field source can be presented in the following 
form: 

1 J'^'^(^) = -87rVi/2 / 5(g«2)<5(gS2)<5[(A;-g)«i]5[(A;-g)Si]e^(*^)+«(<^-''i)) ( + 

(65) 

where the two terms 

^ [jk - g)S.]4"X-' - [{k - q)u,]E[^X^^ - |yr(^2Afe^) 

(k - qniik - q)J:,]' - [{k - q)u,n ^ > 

could be associated with contributions of the first and the second string, respectively. It 

has to be noted, however, that in the second order the axion field is generated collectively 

by a source term containing symmetrically the contributions of the first order coming from 
both strings. It can be checked that the antisymmetric tensors Q^*^ satisfy the transversality 
conditions 

k^^QT = A;^gr = (68) 

ensuring in turn the transversality of the current iJ^'^{k). 

Integration over g-space is performed as follows. First, using 5[{k — q)ui\, we integrate 
over which fixes the value 

= e = uj. (69) 

Then, using the delta-functions 5{qu2) and 5(gE2), one integrates over (f and q^ obtaining 

9 a;cosQ; , a; sin a , 
q = , = ■ 70 

V V 

The remaining delta-function 6[{k — q)^i] does not further fix the vector g^, but imposes an 
extra relation on the parameters : 

k' = (71) 

V 

which is nothing but the Cerenkov condition for radiation emitted by an effective source 
moving with the velocity Vp along the axis. Clearly, this motion is the motion of an 
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effective source localized in the region surrounding the point of minimal separation between 
the strings where their deformation described by the first order interaction is maximal. This 
condition can be rewritten as the usual Cerenkov's cone condition on the angle of emission 

cose = l/vp, (72) 

which makes sense if fp > 1. So the axion radiation is emitted inside the Cerenkov cone 
oriented along the axis (Fig. 2). In the limiting case Vp = 1, radiation is emitted strictly 
along the direction; this occurs when the string inclination angle a is related to the string 
Lorentz factor as 

cosa = 7~\ with 7 = (1 — (uiU2)^^)2 . (73) 

For cos a < 7"^ no radiation is emitted at all. Therefore, string axion bremsstrahlung has 
essentially Cerenkov nature. But, rather surprisingly, it has another interpretation as a 
point charge bremsstrahlung in 2 + 1 electrodynamics. We will come back to this in Sec. 7. 

Integration over is performed using the contour integration in the complex plane. 
Denoting the impact parameter d = d2 — di we obtain for the first term or (j65|) 

-''I'VMq') Trh{-iK,)e--^ ^^^^ 



where /i stands for the remaining integrand, and 



Aci = — . (75) 

7f 



Similarly, the second term in ()65|) gives the integrals 



(76) 



where 

K2 = — , c. = cos a H V. (77) 

V UJ 

Therefore, the whole integral over d'^q of the first term in (j65|) is proportional to the value 
of the integrand in the complex point q^ = q^ where 

(t = — [f, — — cosa, sina], (78) 

with the corresponding value of Xf being 

Xf = ^ [0, icosa, 1/7, 0]. (79) 
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For the second term one finds instead the following quantities 

k"- f/^ = ^[0, z, 1,0], (80) 

X2 = ^ [f , — i cosa, — cos a, sina]. (81) 

Finally, projecting the current on the polarization tensor fl2f)|l . we obtain the following 
lowest order approximation for the quantity (jHU: 

1 J(k) = l^^?!^M5(cos^ - 1/v,) ■ (82) 
uj^'-fv sm & 

(cos — ^7 COS a sin (p) — [cos (p — iiv^^"^ sin 6 + sin (p)] ] , 

where ^ = cos a + f sin 6' sin y?. Two terms here may be attributed to the contributions 
of the first (remaining at rest) and the second (moving) strings; in what follows we set 
/i = /2 = /, /Wi = At2 = /W. The total amplitude is therefore a superposition of two terms 
exhibiting different frequency cutoff. The first term's exponential factor leads to the isotropic 
condition 

^ < (83) 
while the second gives a 9?— dependent condition 

- < ^ = '—r—. (84) 



d{cos a + u sin ^ sin (p) 

In other words, the distribution of radiation on the Cerenkov cone is anisotropic. This 
feature becomes especially pronounced in the ultrarelativistic case. Using the identity 

cos^ a = v^ sin^ 6 + 7"^ (85) 

which holds on the radiation cone, it is easy to show that when 7 — > 00, the quantity ^ has 
a sharp minimum at p = —tt/2 corresponding to the direction of the moving string in the 
rest frame of the first string (see Fig. 2): 

e ~ ^ (1 + /3Vcosa2)], (86) 

27"^ cos a 

where P = tt/2 + p 1. Due to the factor in the second term in the amplitude (jSH), 
radiation is peaked in the direction p = —7t/2 within the narrow region 

/? < 7"'- (87) 
15 



k3 










'/ 

// 







FIG. 2: Radiation peaking on Cerenkov cone in the direction of motion of the relativistic sting. 
The exponential factor exhibits similar peaking, provided 

7cosa>l, (88) 

so the frequency range extends up to 



~ d cos a 



in the angular region (jSTj). 
At the Cerenkov threshold 

V = sin a, 

the cone of emission shrinks 6^0, in this case 

= cos a = 



(89) 



(90) 



(91) 



and the amplitude ()82|1 can be considerably simplified. Since in this limit kdn —>■ 0, the 
oscillating factors exp (ikdn) become equal to unity, and one can check that the leading 
contributions of two terms in (jH^ cancel while ^ 0, that is the whole expression does not 
diverge. The remaining finite contribution reads 



,J(k) 



-e— '^M5(cose- 1) 



(92) 



Inserting this into the Eq. (j3Uj) and introducing the normalization length / via the relation 



^^(cos^- 1) = — (5(cos0- 1) 
2% 



(93) 
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one obtains after integration over angles in d^k = uo"^ duo d cos 9 dip the total radiated energy 
per unit length of the string at rest 

- ' ^ ' ■exp(-2 — )duj, (94) 



where A is the frequency cutoff introduced in order to remove the infrared divergence of the 
spectral distribution of radiation. At the high frequency end the spectrum extends till the 
boundary ()8H|) . Integration over frequencies is performed in terms of the integral exponential 
function 

For d <^ vyA, 

where C is the Euler constant, e*" = 1.781072418. For d ^ vfA, one has 

f 2d \ 2d f 2d \ 

in this case radiation is exponentially small. A particular value, Ei(l, 1) = 0.219, can be 
used to estimate the energy loss for intermediate impact parameters. For d = vfA one has 



— = 2148 ^„ 
I fj,^ 



— = 2148^. (98) 



For generic values of parameters a, f it is convenient to rewrite the amplitude in terms 
of the quantity k = 7 cos a which is invariant under world-sheets reparameterizations char- 
acteristic of the relative motion of the strings. Setting without lack of generality di = we 
obtain: 



iJ(k) = <^exp (cosy^-msiny?) 

cos if — in (^K sin y? + ^/ k"^ — 1^ 



— exp 



ud 



(99) 



(k + \/ ~ 1 sin (pj 

This is still too complicated to perform the subsequent integration analytically, so we pass 
to the limit k — > 00, assuming f — > 1 and the inclination angle a not close to the threshold 
of shrinking cone. Then the second term exhibits peaking around the direction (p = —tt/2 
which marks the plane of motion of the relativistic second string (recall that the first string 
is at rest in the Lorentz frame chosen). Keeping only the second term which is dominant in 
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FIG. 3: The spectral function F{z). 

the main frequency range we obtain the spectral-angular distribution in the vicinity of this 
direction as follows 



1 dP^ QAtt^Pk^I- K^0 



The spectrum exhibits an infrared divergence, and in the forward direction = it extends 
up to a; ~ Wmax, where 

t^max = (101) 

Choosing as an infrared cutoff an inverse length parameter A and integrating and over 
frequencies we obtain the angular distribution of the total radiation 



= T. 7-. — -— tIIjI 1, -. 



/ d(5 



:i + fi:2/3 



7kA 



(102) 



Since the integral exponential function decays exponentially for large argument, the total 
radiation is peaked around /5 = within the angle 

13 (103) 

One can also obtain the spectral distribution of radiation by integrating p()()|) over the 
variable (3, which can be extended to the full axis in view of the exponential decay of the 
integrand. The result reads: 



IdP^ IQiT^fd^, , uod 
= —F{^)^ z = — 



I duj 



37/^2 



where 



F = (8z^ + 8z + 6)^= + (8^^ + 122 + 6 - 3/z)(erf(v^) - 1) 



(104) 



(105) 
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FIG. 4: Dependence of radiation energy on the impact parameter, y = d/{'~fKA). 

The spectrum is shown on the Fig. 3. For small frequencies u <C 'yn/d the function F{z) 
has a logarithmic divergence 

F{z) ~ -, (106) 



while for large z it decays exponentially 

F(z) 



12e- 



(107) 



Finally, the total radiation per unit length is obtained by integrating ()104|) over frequencies 
using a cutoff length A: 

(108) 



—(I){y), y 



where 



M = 12 



3/i2 

113 3 



7r'^H2'2'2'2' ^ 



7kA ■ 



3 In {Aye^) + 



7 



(109) 



where the generalized hypergeometric function is used. The function (j){y) is plotted on Fig. 
4, for small y the leading term is the logarithm, while the asymptotic value for large y is 
0(oo) = 7/2. 

The speed of the effective radiation source Vp tends to infinity when the two strings 
become parallel. In this case (corresponding to k = 7 in the above formula) the Cerenkov 
cone opens up to 6 = tt/2, and the whole picture becomes essentially 2 + 1 dimensional. 
Since this case is in fact generic (superluminal strings can be always "parallelized" by the 
world-sheet reparameterizations ()43|1 and suitable space-time Lorentz transformations), this 
gives an alternative description of the effect as bremsstrahlung under collision of charges in 
2 + 1 electrodynamics. 
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VI. BREMSSTRAHLUNG IN 2+1 ELECTRODYNAMICS 



Classical system of parallel strings interacting via the antisymmetric two-form potential 
B^i, can be equivalently presented as the system of point charges interacting via the Maxwell 
field in 2+1 dimensions (in what follows in this section the Greek indices run over 
fi = 0,1,2). Indeed, substituting the constraint into the action Q and assuming that all 
variables are a;^-independent, one can rewrite the action as 

S = -J2[ ("'-^^^ + ^n^M^O - \ [ F^^F'^-'d'x. (110) 

The correspondence between the quantities here and in the previous section looks as follows: 



^^ = ^2/53^, m„ = /i„/, e„ = 2V2/7r/„, (111) 

and F^y = d^Ay — d^A^, where / is the normalization length along the suppressed coordinate 
x^. The equations of motion in the gauge d^A^ = read 

^n^n ~ ^nF^ ui-ni (H2) 

□A^ = -J/^, (113) 

where = + J!^, 

Jii = en <5^(X - Xn{Tn))dTn, (114) 



and only mutual interaction should be taken in the particle equations of motion after the 
classical mass renormalization. 

Performing standard calculations we obtain for the energy loss under particle collision: 

1 /■ . „ 



/ '''^\^ik')\J,{k)\'d'k. (115) 



In 2+1 dimensions a photon has the unique polarization state, so alternatively this expression 
can be rewritten (after integration over k^) as 

1 r kf^ 

P' = ^2 J j^\Jm'd\ (116) 

where we put = |k|, = (0,e), J(k) = J^{k)e^, and two-dimensional vectors are 
parameterized as 

k = u; [cos (p, sin ip], 

e = [— siny?, COS9?]. (117) 
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Perturbation theory with respect to the formal charge parameter e (whose dimension in 
the 2 + 1 electrodynamics is Z"^/^ in the units h — c— 1) is constructed in the same way as 
in the Sec. II, introducing the input world-lines 

o< = < + <r„, (118) 

with space-like impact vectors and unit time- like 3- velocities u^: 

< = [1,0,0], u^,='y[l,0,-v]. (119) 

After the first iteration one obtains 

where d = d2 — di, and similarly for the second charge. The first order current, which gives 
the lowest order contribution to radiation, reads 



J^^Y.^- f - < l<^-) ^ - 0^r^{rn)) • (121) 

n=l,2 



Its Fourier-transform is 

6162 
2n 



f 6{qu2)6[{k - g)Mi]e*('^i^+«'^) (^Qt + ^Qi^) d^q, (122) 
J \mi 1712 J 



where 



_ {kuifui^ - {kui){ ku2)u1 + {uiU2)[{kq)u'l - {kui)q>'] 



2 



q^{kui) 

_ {ku2Yu'i - {kui){ku2)u^ + {uiU2)[{kq)u^ - {kU2){k - qf] 

{k-qnku2f ■ ^^^^^ 

Integration over q is performed similarly as in the case of strings: two delta-functions are 
used to integrate over g^, the remaining integral over q^ is evaluated using a contour 
integration in the complex plane. Finally, projecting the amplitude on the polarization 
vector e'' we obtain 

6162 / Q^Q-'^d/vj+ikdi ^^Q-uid^/v+ikd2 \ 

iJ(k) = {-cosip + i'ysiiKp) — — [cos (p - i{v^'y'^ + sin ip)] , 

2u!^v \ mi m27^ / 

(124) 

where ^ — 1 + vsin(f. 

As we argued in the Sec. 3, our iteration scheme converges if the relative velocity is 
relativistic, and for cosmological applications one is interested by moderately relativistic 
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velocities when 7 is not large. However, for general 7 it is difficult to perform computation 
of the total energy loss analytically, so we continue calculations in the ultrarelativistic case 
7^1 when additional simplifications can be made. Hopefully this should give a reasonable 
approximation for 7 of the order of several units. One can show that then the second term in 
(I124j) is dominant (the asymmetry of contributions from the first and the second charges is 
due to our choice of the reference frame: the first charge is at rest before the collision, while 
the second is moving with the relativistic velocity). The spectral distribution of radiation 
due to the second term extends till lu ~ in the narrow angular region 6(f ^ I/7 around 

the direction ip = — 7r/2, while the first term has a lower (~ "y/d) frequency cutoff. The total 
contribution of the second term to the radiation loss is greater than that of the ffist one by a 
factor of 7, so in what follows we neglect the first term. Setting ei = 62 = e, mi = m2 = m 
we find from the Eq. ()116|) the spectral-angular distribution of radiation as follows 

dP° _ 1 [cos^ ip + (sin ^ + ^^y] e'^'^'^ 

dud^ ~ 32^^ ^^76^4 • ^ ) 

For 7^1, radiation is peaked in the direction of motion of the relativistic charge = — 7r/2, 
and in the vicinity of this direction (cf. ()86|l ) 

+ (126) 



where /3 -C 1, and 



exp . (127) 



dujdf3 SttWcu (1 + 72/32)4 ^2 
The spectrum has an infrared divergence, and in the forward direction /3 = it extends up 



to ~ cUmax, where 



72 



Integrating over frequencies from the infrared cutoff A ^ we obtain 



d(3 87r2m2 (1 + ^2^2)4 \^ ' 

Taking into account relations (jlllj) between parameters of 3+1 and 2+1 problems, one can 
see that this expression coincides with ()102j) for 7 = k, i.e. for a = 0: 

1 dP' 647r VV (1 - L d^ + f^\ 



/ d[3 /i2 (1+^2^2)4 \^ ' 
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VII. COSMOLOGICAL ESTIMATES 



In the standard axion cosmology scenario, which assumes the Peccei-Quinn phase transi- 
tion after inflation, the global string network is formed at the temperature Tq ~ /. Strings 
initially move with substantial friction [l^, [3, ll^ due to scattering of particles present 
in the hot cosmic plasma. At some temperature T* < / the string network enters the 



scale-invariant regime |20|, |21L |22| when strings move almost freely with relativistic veloci- 
ties 23|, |2^ and are able to form closed loops via interconnections. Finally, at much lower 
temperature Ti of the QCD phase transition, axions become massive and strings disappear. 
Our mechanism presumably should work in the temperature interval (X,,, Ti), though an 
additional contribution from the damped epoch Tq < T < can also be expected. The 
scaling density of strings is determined from numerical experiments, it can be presented as 



Ps = (131) 



with ( varying from 1 to 13 in different simulations. 

Our aim is to estimate the total energy density Ea of axions produced via the 
bremsstrahlung mechanism during the cosmological period from t^, (corresponding to tem- 
perature T^) till the QCD phase transition moment ti. These are commonly estimated 
as Q 

t.(sec) ~ 10-^° (TTdT-Tj) ' (132) 



^IQi^GeV^ 

Consider the scattering of an ensemble of randomly oriented straight strings on a se- 
lected target string in the rest frame of the latter. Since the dependence of the string 
bremsstrahlung on the inclination angle a is smooth, we can use for a rough estimate the 
particular result obtained for the parallel strings (a = 0) introducing an effective fraction u 
of "almost" parallel strings (roughly 1/3). Then, if there are strings in the normalization 
cube = L^, we have to integrate the radiation energy released in the collision with the 
impact parameter d = x over the plane perpendicular to the target string using the measure 
N / L"^ ■ 2TTxdx. Actually we need the radiation power per unit time, so, for an estimate, we 
have to divide the integrand on the impact parameter. This quantity should be multiplied 
by the total number of strings to get the radiation energy released per unit time within 
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FIG. 5: Dependence of the axion radiation rate on the Lorentz factor of the colhsion ^{w), w 
L/(7^A). 



the normahzation volume. Therefore, for the axion energy density generated per unit time 
we obtain: 

As an expression for P° we will use here the Eq. p08|) . and we have to integrate over the 
impact parameter. Taking into account that the string number density is related to the 
energy density ()131|) via 



V 



(135) 



introducing the integral 



aw) 



ct){y)dy = 12J- 



'0 

— 2F2 



2-^2 



1 1.11. 

'2'~2' 2' 2' 



-w 



1 1.13. 

2'~2' 2' 2' 
13 



-w 



3 In [Awe^] w, 



(136) 



and using for /i the expression (|TT|) with only the second (leading) term, we obtain 



dfa 
dt 



where 



w 



72 A' 



(137) 



(138) 



Since in the cosmological context L A t, one can take w = 7^^ for an estimate. 

The realistic value of 7 is of the order of unity, while our formulas were obtained in the 
7^1 approximation, since we were keeping only the second term in the exact amplitude 
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()99|) . But an independent calculation shows that for 7 ~ 1 the contribution from the first 
term is of the same order, so we can hope to get reasonable order of magnitude estimate 
using the above formulas for 7 not large. The full 7-dependence is given by the function 
7^^(7~^), where ^{w) is plotted on Fig. 5. For 7 = 

I'ar') - 8.137, (139) 

while for 7 = 5 

73^(7-2) ~ 55.83. (140) 

Another source of uncertainty is the numerical value of the parameter (, which was obtained 
in different simulations within the region between 1 and 13. Now we have to integrate the 
radiation power over time from to ti as given by ()132|1 and ()133|) . Using the integral (with 
substitution x = 1/t) 

2 

^ dx=-:^+2U{x^), (141) 



In^ X In X 

where the logarithmic integral function is introduced, which for large x is equal approxi- 
mately to 

we find the energy density of axions at the moment ti ^ t*: 

~ ^ 2(. 1. V (143) 



where the ratio t^f/to can be estimated as 

2 

(144) 



Finally, dividing by the critical energy density at t = ti 



3'm2 

^cr = 4^, mpi = 1.22 ■ lO^^GeV, (145) 
we obtain the following estimate for the relative contribution of bremsstrahlung axions 

f!?~0.5.10-(i^)''(^)"'\ (146) 

where we have used the value (jl39|) . This estimate shows that our new mechanism is very 
sensitive to the value of / and gives an upper bound on axion window of the order of 
lO^^GeV. We postpone more precise cosmological estimates for a separate publication. In 
particular, more careful analysis is needed to find radiation loss for moderately relativistic 
string velocities. 
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VIII. CONCLUSIONS 



In this paper, we have suggested a new mechanism for the axion emission by the global 
string network: the bremsstrahlung under string collisions. As far as we are aware, this 
effect was not discussed in the context of the cosmic string scenario so far. Though it is of 
the second order in the axion coupling constant, rough cosmological estimates show that it 
is not small and gives a contribution of the same order of magnitude as radiation due to 
oscillating loops. 

We have found the radiation amplitude using the perturbation theory whose validity 
is restricted by relativistic (though not necessarily ultrarelativistic) velocities of colliding 
strings. The frequency spectrum has an infrared divergence which is not an artifact of our 
approximation, but rather is the effect of the space-time dimensionality: as we have shown, 
an equivalent description of the axion bremsstrahlung from strings in 3 + 1 dimensions is 
provided by the electrodynamics of point charges in 2 + 1 dimensions, where its origin lies 
in the logarithmic dependence on distance of the Coulomb potential. Unfortunately the 
integration over the spectrum and the angular distribution can be performed analytically 
only in the ultrarelativistic case, and the final formulas obtained relate to this situation. 
But, in view of a smooth dependence of the radiation efficiency on the Lorentz factor of the 
collision, we hope that our results give a reasonable estimate of the effect for values of the 
Lorentz factor of the order of several units as well. 

Prom a purely theoretical viewpoint two aspects are worth to be discussed. The first is 
the Cerenkov nature of radiation which is associated with the possibility of the superluminal 
motion of an effective source located around the point of minimal separation between the 
strings. In perturbative setting, the strings, which are straight in zero order approximation, 
get deformed under axion interaction, with the deformation being maximal near this point. 
This region therefore acts as an effective source of radiation arising in the next order ap- 
proximation. This explains why radiation is concentrated on the Cerenkov cone directed 
along the trajectory of the effective source. Another interesting feature is that the effect has 
an alternative 2-1-1 interpretation as the bremsstrahlung of point charges. This is due to 
the fact that two strings inclined with respect to each other and moving in the superluminal 
regime can be made parallel by suitable reparamctcrizations of their world-sheets and the 
choice of the space-time Lorentz frame. This transformation is only feasible once the relative 
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string motion is superluminal in the sense described above. Classical dynamics of parallel 
strings interacting via the axion field in 3 + 1 dimensions is equivalent to dynamics of point 
charges in 2 + 1 electrodynamics, therefore the string axion bremsstrahlung is reduced to 
the point particle bremsstrahlung in lower dimension. 
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